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^ : Abstract 

It is shown that the expected amplitudes and specific correlation properties of 

■ the relic (squeezed) gravitational wave background may allow the registration 
of the relic gravitational waves by the first generation of sensitive gravity- wave 
detectors. 
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Relic gravitational waves are inevitably generated by strong variable gravitational field 
of the very early Universe which parametrically (superadiabatically) amplifies the zero-point 
quantum oscillations of the gravitational waves [1] . The initial vacuum quantum state of each 
pair of waves with oppositely directed momenta evolves into a highly correlated multiparticle 
state known as the two-mode squeezed vacuum quantum state [2]. (For a recent review of 
squeezed states see [3]). In cosmological context, the squeezing manifests itself in a specific 
standing- wave pattern and periodic correlation functions of the generated field [4]. The main 
point of this paper is to show that this signature may significantly facilitate the detection of 
the relic (squeezed) gravitational wave background. It is possible that the appropriate data 
processing may allow the detection of relic gravitational waves by the forthcoming sensitive 
instruments, such as the initial laser interferometers in LIGO [5], VIRGO [6], GEO600 [7]. 

Before analyzing the theoretical predictions in light of the current experimental situation, 
it is necessary to explain the origin of these predictions. 

We consider the cosmological gravitational wave field defined by the expression 

ds 2 = a 2 {r])[dr] 2 - (5 tj + h^dx'dx 1 ] . (1) 
The Heisenberg operator for the quantized real field hij can be written as 
(' f - . ^ ••>• , , 1 

s=l 



^' x) = (2^7* L/ n ^ PiM) ^m i hMeim Cn + ^ )e ~ mx c «J ' (2) 



where C = VTQtt I pi and the Planck length is Ipi = (Gh/c 3 ) 1 ^ 2 , the creation and annihilation 
operators satisfy c n |0) = 0, [c n , c m ] = 5 ss >5 3 (n — m), and the wave number n is related to 
the wave vector n by n = (dijii 1 ^) 1 ^ 2 . The two polarisation tensors Pij(n) (s = 1,2) obey 
the conditions 



s 



Pijn 1 = 0, P i:j S lJ = 0, P i:j P 13 = 25 ss ,, Pij(-n) = Pijin). 

s 

For every wave number n and each polarisation component s, the functions h n (v) have the 
form 

L(v) = ^r[un(v)+Z(v)], (3) 
a{rj) 

where u n {rj) and v n {rj) are expressed in terms of the three real functions: r n - squeeze 
parameter, n - squeeze angle, 9 n - rotation angle, 

u n = e i9n cosh r n , v n = e ~ i{ - 6n ~ 24>n) sinh r n . (4) 

The functions r n {rj), <f> n (v)i 9n(v) are governed by the dynamical equations [4]: 

r' n = — cos20 n , <p' n = —n sin20„coth 2r n , 9' n = —n sin20 n tanh r n , (5) 

(X (X (X 

where ' = d/drj, and the evolution begins from r n = 0, which characterizes the initial vac- 
uum state. The dynamical equations and their solutions are identical for both polarisation 
components s. 
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The present day values of r n and 4> n are essentially all we need to calculate. The mean 
number of particles in a two-mode squeezed state is 2sinh 2 r„ (for each s). This number 
determines the mean square amplitude of the gravitational wave field. The time behaviour 
of the squeeze angle (p n determines the time dependence of the correlation functions of the 
field. The amplification (that is, the growth of r n ) governed by (5) is different for different 
wave numbers n. Therefore, the present day results depend on the present day frequency v 
[y = cn/2na) measured in Hz. 

In the short- wavelength (high-frequency) regime, that is, during intervals of time when 
the wavelength \(n) = 2-na/n is shorter than the Hubble radius l(n) = a 2 /a', the term n 
in (5) is dominant, and the functions (f) n (r)), n (n) are 0„ = — n(n + const), 6 n = <p n . The 
factor cos 2<f> n is a quickly oscillating function of time, so the squeeze parameter r n stays 
practically constant. 

In the opposite, long-wavelength regime, the term n can be neglected. The function <p n 
is ta.n (f) n (rj) ~ const/ a 2 (n), and the squeeze angle quickly approaches one of the two values 
4> n = or 4> n = 7r (an analog of the "phase bifurcation" [8]). The squeeze parameter r n (n) 
grows with time according to 

r n (n)^ln^l , (6) 

where a* is the value of a(n) when the long-wavelength regime, for a given n, begins. The 
final amount of r n is 

r n ps In— , (7) 

where a** is the value of a(n) when the long-wavelength regime and amplification come to 
the end. 

After the end of amplification, the accumulated (and typically large) squeeze parameter 

S S * 

r n stays approximately constant. The complex functions u n (n) + v n (n) become practically 

s s * 

real, and one has h n (v) ~ h n (v) ~ ^e r ™ cos 0„(?7). Every mode n of the field (2) takes the 
form of a product of a function of time and a (random, operator-valued) function of spatial 
coordinates, that is, the mode acquires a standing-wave pattern. The periodic dependence 
cos 4> n (n) is a subject of our further inquiry. 

The numerical results depend on the concrete behaviour of the pump field represented 
by the cosmological scale factor a(n). We know that the present matter-dominated stage 
a(n) oc rf was preceeded by the radiation-dominated stage a(n) oc n. The function a(n) 
describing the initial stage of expansion of the very early Universe (before the era of pri- 
mordial nucleosynthesis) is not known. It is convenient to parameterize a (n) by power-law 
functions in terms of n, since they produce gravitational waves with power-law spectra in 
terms of v [1]. Concretely, we take cii(n) at the initial stage of expansion as aiirf) = l \i]\ 1+lS 
where n grows from — oo, and /3 < — 1. From n — rji, r\\ < 0, the initial stage is followed by 
the radiation-dominated stage a e (n) = l cie(v — Ve) and then, from n = i]2, by the matter- 
dominated stage a m (n) = lo^mijI^Vm) 2 ■ The constants a e , a m , n e , r\ m are expressed in terms 
of the fundamental parameters l Q , f3 through the continuous joining of a(n) and a'(n) at rji, 
1]2- The present era is defined by the observationally known value of the Hubble radius 
In = c/H ps 2 x 10 28 cm. We denote this time by n R and choose n R — r\ m = 1, so that 
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a (VR) = 21 h- The ratio a(r]R) / '0(772) = z is believed to be around z = 10 4 . Some information 
about l and /3 is provided by the data on the microwave background anisotropies [9, 10], 
and we will use it below. 

The known function a(rj) and equations (5) allow us to find the present day values of 
r n as a function of n. Let the wave numbers nninn = 47r), n m {n m = y/znn), n c denote 
the waves which are leaving the long-wavelength regime at, correspondingly, t)r, 772, 771- 
The n c is the value of n for which a** = a*. The shorter waves, with n > n c , have never 
been in the amplifying long- wavelength regime. (The present day frequency v c is around 
10 10 Hz, as we will see below). Thus, r n = for n c < n, r n = /n[(n/n c ) /3 ] for n m < n < n c , 
r n = ln[(n/n c ) /3 ~ 1 (n m /n c )] for n H < n < n m , and r n = ln[(n/n^ +1 (n m /n H )(n c /n H )\ for 
n < n-H- The e Tn is much larger than 1 for all frequencies n < n c < u c ). 

The mean value of the field hij is zero, (0|/iy(?7,x)|0) = 0. The variance 

(0|/iyfax)/i«(i7,x)|0) = (/i 2 ) 

is not zero, and it determines the mean square amplitude of the generated field - the quantity 
of interest for the experiment. Taking the product of two expressions (2), one can show that 



£*2 »oo 2 g 2 £*2 »oo 

(h 2 ) = — nJ2\hn(v) dn = ^T^2 ndn(cosh2r n + cos20 n sinh2r„). (8) 

ZiTT JO -1 7T Cl JQ 



C 2 r°° A \s . . 2 , c 2 • - 

s=l 

Eq. (8) can also be written as 

(h 2 )= fVfn,^, 
Jo n 

where, for the present era, 



C 1 

h(n,i]) ~ — -ne Tn cos <f) n {i]) . (9) 

7T a(r) R ) 



The further reduction of this formula gives 

/^X/3+2 



h(n) « A — , n<n H , (10) 

/ n \P 

h(n,rj) ps A cos n (77) — , n H < n < n m , (11) 
/ n {ri}{\ 

h(n, i]) ^ A cos 0„(?7) — — , n m < n < n c , (12) 

Vn H / Vn m / 

where 

4 80F2^+ 2 



The available information on the microwave background anisotropies [9, 10] allows us to 
determine the parameters A and /3. The quadrupole anisotropy produced by the spectrum 
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(10) - (12) is mainly accounted for by the wave numbers near %. Thus, the numerical value 
of the quadrupole anisotropy produced by relic gravitational waves is approximately equal 
to A. Since (according to [11]) the quadrupole contribution of relic gravitational waves is 
not smaller than that produced by primordial density perturbations, this gives us A m 1CT 5 . 
We do not know experimentally whether a significant part of the quadrupole signal is indeed 
provided by relic gravitational waves, but we can at least assume this. The evaluation of 
the spectral index n of the primordial perturbations resulted in n = 1.2 ± 0.3 [10] or even 
in a significantly higher value n = 1.84 ± 0.29 [12] (see also [19], where one of the best 
fits corresponds to n = 1.4). We interpret [13] these evaluations as an indication that the 
true value of n lies somewhere in the interval n = 1.2 ~ 1.4 (hopefully, the planned new 
observational missions will determine this index more accurately). Since n = 2/5+5, this gives 
us the parameter (3 in the interval /3 = —1.9 ~ —1.8. In fact, the value n = 1.4 (j3 = —1.8) is 
the largest one for which the entire theoretical approch is well posed, since this value requires 
the pumping field to be too strong: the Hubble radius at the end of the initial stage would 
have been only a little larger than l P i. With the adopted A and (3, the frequency v c falls 
in the region v c = 10 10 Hz or higher. An allowed intermidiate stage of expansion governed 
by a "stiff" matter [14] can affect the spectral slope at frequencies somewhat lower than u c , 
but still outside the interval accessible to the existing detection techniques. The details of 
a short transition from the initial stage to the radiation-dominated stage are irrelevant for 
our discussion, since they can alter the signal only at frequencies around v c . It is worth 
recalling [13] that a confirmation of any n > 1 {J3 > —2) would mean that the very early 
Universe was not driven by a scalar field - the cornerstone of inflationary considerations - 
because the n > 1 (f3 > —2) requires the effective equation of state at the initial stage to 
be e +p < 0, but this cannot be accomodated by any scalar field with whichever scalar field 
potential. Obiously, the available data do not prove yet that n > 1 but this possibility is not 
ruled out either. It is important to emphasize that the existing disagreement (see [11]) with 
regard to the validity of the inflationary prediction of infinitely large amplitudes of density 
perturbations for the part of the spectrum with the spectral index n = 1 (f3 — —2) does not 
affect the assumptions and conclusions of this paper which analyses the cases (3 > — 2. 

We switch now from cosmology to experimental predictions in terms of laboratory fre- 
quencies v and intervals of time t (cdt = a(i] R )dri). Formula (12) translates into 



where Vh — lO~ 18 Hz, and t v is a deterministic (not random) function of frequency which 
does not vary significantly on the intervals Av ~ v. We take v = 10 2 Hz as the representative 
frequency for the ground-based laser interferometers. The expected sensitivity of the initial 
instruments at v — 10 2 Hz is h ex = 10~ 21 or better. The theoretical prediction at this 
frequency, following from (13), is h t h = 10~ 23 for f3 = —1.8, and h t h = 10~ 25 for (3 = —1.9. 
Therefore, the gap between the signal and noise levels is from 2 to 4 orders of magnitude. 
This gap should be covered by a sufficiently long observation time r. The duration r depends 
on whether the signal has any temporal signature known in advance, or not. 

It appears that the periodic structure (13) should survive in the instrumental window of 
sensitivity from v\ (minimal frequency) to u 2 (maximal frequency). The mean square value 
of the field in this window is 




(13) 
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/ tffa t )— = 10~ 14 —— / cos 2 [27rz/(t - t,)]z/ 2/3+1 dz/ . (14) 

Because of the strong dependence of the integrand on frequency, z/~ 2 - 6 or z/ -2 - 8 , the integral 
(14) is determined by its lower limit. This gives 

rV2 Ail / \ 2/3+2 

/ h 2 (lS,t)— « 10- 14 ( — ) COS 2 [27TI/i(t-ti)] . (15) 

The explicit time dependence of the variance of the field, or, in other words, the explicit 
time dependence of the (zero-lag) temporal correlation function of the field, demonstrates 
that we are dealing with a non- stationary process (a consequence of squeezing and severe 
reduction of the phase uncertainty). Apparently, the search through the data should be 
based on the periodic structure at v — v\. 

The response of an instrument to the incoming radiation is s(t) = Fijh^ where Fij 
depends on the position and orientation of the instrument. The cross correlation of responses 
from two instruments (0|si (t)«2 (t) |0) will involve the overlap reduction function [15 - 18], 
which we assume to be not much smaller than 1 [17]. The essential part of the cross 
correlation will be determined by an expression of the same form as (15). 

The signal to noise ratio S/N in the measurement of the amplitude of a signal with 
no specific known features increases as (rz/i) 1 / 4 . If the signal has known features exploited 
by the matched filtering technique, the S/N increases as (rz/i) 1 / 2 . The guaranteed law 
(tvi) 1 / 4 requires a reasonably short time r = 10 6 sec to improve the S/N by two orders of 
magnitude and to reach the level of the predicted signal with the extreme spectral index 
(3 = —1.8. If the law (r^) 1 / 2 can be implemented, the same observation time r = 10 6 sec 
will allow the registration of the signal with the conservative spectral index (3 = —1.9. Even 
an intermediate law between (rz/i) 1 / 4 and (rz/i) 1 / 2 may turn out to be sufficient. For the 
network of ground-based interferometers the expected v\ is around 30Hz, but we have used 
V\ = 10 2 Hz for a conservative estimate of r. If the matched filtering technique can indeed 
be used, it can prove sufficient to have data from a single interferometer. 

For the frequency intervals covered by space intereferometers, solid-state detectors, and 
electromagnetic detectors, the expected results follow from the same formula (13) and have 
been briefly discussed elsewhere [13]. 

In conclusion, the detection of relic (squeezed) gravitational waves may be awaiting only 
the first generation of sensitive instruments and an appropriate data processing strategy. 

I appreciate useful discussions with S. Dhurandhar and B. Sathyaprakash. 
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